Energy transfers and locality in magnetohydrodynamic turbulence 
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The shell-to-shell energy transfer rates for magnetohydrodynamic (MHD) turbulence are com- 
puted analytically, which shows local energy transfer rates from velocity to velocity, velocity to 
magnetic, magnetic to velocity, and magnetic to magnetic fields for nonhelical MHD in the inertial 
range. It is also found that for kinetic-energy dominated MHD fluid, there is a preferential shell-to- 
shell energy transfer from kinetic to magnetic energy; the transfer is reversed for magnetic-energy 
dominated MHD fluid. This property is the reason for the asymptotic value of Alfven ratio to be 
close to 0.5. The analytical results are in close agreement with recent numerical results. When mag- 
netic and kinetic helicities are turned on, the helical contributions are opposite to the corresponding 
nonhelical contributions. The helical energy transfers have signiflcant nonlocal components. 

lO ' PACS numbers: 47.65. +a, 52..30.Cv, 52..35.Ra 
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' Turbulent fluid and plasma flows exhibit complex behaviour. One such phenomena is the energy transfers among 
various scales. For fluid turbulence the energy transfer issues have been investigated in great details. However, 
detailed analysis of these processes is lacking in magnetohydrodynamic (MHD) turbulence. Detailed understanding 
■ of energy transfer is useful for understanding various physical process, for example, dynamo mechanism to generate 
magnetic fleld in astrophysical objects. These results are also useful in modelling MHD flows and in simulations. For 
example, we need to model backscatter and forward energy transfer for large-eddy simulations. In the present paper 
we investigate the above issues analytically. 

Kolmogorov's fluid turbulence phenomenology for incompressible turbulence is based on local energy transfer be- 
Ph \ tween wavenumber shells. There are several quantitative theories in fluid turbulence about the amount of energy 
1—1 • transfer between neighbouring wavenumber shells. For examples, Kraichnan Q showed that 35% of the energy flux 

1 ^ i ' comes from wavenumber triads where the smallest wave-number is greater than one-half of the middle wavenumber. 

This phenomenology has been verifled using numerical and analytical methods p, 3, 4, 5]. Debliquy et al. recently 
' studied the issues of energy transfers in decaying magnetohydrodynamic (MHD) turbulence using direct numerical 
K>" , simulation (DNS). Alexakis et al. and Mininni et al. Q performed the similar calculations for forced MHD tur- 
)^ • bulence for both heHcal and nonhelical flows. They found that typically, the shell-to-shell energy transfer is local. In 
the present paper we compute the above quantities analytically, and compare them with the numerical results. 

The interactions in MHD are through (u(k), u(p), u(k — p)) and (b(k), u(p), b(k — p)) triads, where (u,b) are 
the velocity and magnetic fields respectively, and k, p, and k — p are the wavenumbers of the triad. Kraichnan Q 
. gave a general formalism to compute the magnitudes of triad interactions using transfer function S{k\p, q) In this 
' paper we will compute the shell-to-shell energy transfer using a modified method called mode-to-mode energy transfer 
, rate S{k\p\q), which represents the energy transfer mode from p to mode k, with mode q acting as a mediator. The 
• ^ new formalism is necessary for computing the shell-to-shell energy transfer because the earher formaHsm suffers from 
ambiguity arising due to the third leg of the interaction (see Jil, llj). The calculation is done using perturbative 
!• • field-theory up to first-order in perturbation. We take Kolmogorov's spectrum for the energy spectrum as discussed 
!^ ' in current numerical and analytical papers iVi ll4L ll^ iWL Il7l | . Note that the field-theoretic calculations has a lot 

of similarity with Eddy-damped quasi-normal Markovian (EDQNM) approximation. 
$H ' MHD turbulence involves interactions among velocity and magnetic modes, hence energy transfer takes place 
between velocity to velocity, magnetic to magnetic, velocity to magnetic, and magnetic to velocity modes. Debliquy 
et al. 01 computed the shell-to-shell energy transfers in decaying MHD turbulence using simulation data on 512^. 
In their calculation cross helicity, magnetic helicity, and kinetic helicity are negligible. They also took the mean 
magnetic field to zero. Debliquy et al. found forward and local energy transfer from velocity to velocity, and magnetic 
to magnetic fields. Regarding the velocity-to-magnetic energy transfer, for the Alfven ratio greater than approximately 
0.4, the energy transfer is from kinetic to magnetic; the transfer direction is reversed for Alfven ratio less than 0.4. 
Dar et al. Iflj, Alexakis et al. J] and Mininni et al. Q have done the similar analysis for 2D and 3D forced MHD 
turbulence, with forcing at small wavenumber of velocity. They find local energy transfer for velocity to velocity 
fields, and magnetic to magnetic fields in the inertial range. However, the small-wavenumber velocity shells provide 
energy to the small-wavenumber magnetic shells, as well as to the inertial range magnetic shells. We will show in this 
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paper that our theoretical results on inertial range energy transfers are in general agreement with the above numerical 
results. 

Pouquet et al. were the first to investigate whether interactions in MHD turbulence are local or not. Their 
analysis is based on eddy-damped quasi-normal Markovian (EDQNM) calculation. They claimed that nonlocal inter- 
actions exist in MHD due to the mean magnetic field (Alfven effect) and helicity. According to Pouquet et al., the 
local interactions cause the energy cascade, but the nonlocal ones lead to an equipartition of kinetic and magnetic 
energy. In the present paper we will also show that helicity induces nonlocal energy transfers. 

A detailed picture of energy transfers is very useful for understanding turbulence and its modelling. In this paper 
we will show how we can use our theoretical results to argue why the asymptotic state of MHD turbulence is close 
to 0.5. The detailed shell-to-shell energy transfer also provide us important ideas for large-eddy simulations (see 
Debliquy et al. @ for connection with large-eddy simulations) and EDQNM calculation, which assumes local energy 
transfers among wavenumber shells. 

It is well known that compressible turbulence involves energy transfers from pressure fiuctuations to the velocity and 
magnetic field [19 , 21, J2, _2^J . The theoretical, numerical, and observational studies show that the energy spectrum 
deviates from Kolmogorov's spectrum. For example. Burgers equation, which represents fully compressible fiuid, has 
energy spectrum proportional to k~'^. The theory of compressible turbulence is not yet developed as much as that for 
incompressible turbulence. Due to the uncertainty of energy spectrum and other properties in the inertial range, in 
this paper we have confined ourselves to the study of shell-to-shell energy transfer for incompressible turbulence only. 

The outline of the paper is as follows: In Sec. 2 we compute the shell-to-shell energy transfer rates for nonhelical 
and helical MHD. In Sec. 3 we use our results to show why the asymptotic state of MHD turbulence has Alfven ratio 
close to 0.5. The last section. Sec. 4, contains conclusions. 



II. CALCULATION OF THE SHELL-TO-SHELL ENERGY TRANSFERS 

In MHD turbulence, velocity (u) and magnetic fields b interact with each other and among itself to produce complex 
energy transfers. The energy exchange can take place between a u Fourier mode to u Fourier mode, between a b 
Fourier mode to b Fourier mode, or between a u Fourier mode to b Fourier mode. These transfers are studied using 
Kraichnan's formula S'(k|p,q) [9] or mode-to-mode energy transfer rate S'(k|p,q). In this paper we use the mode-to- 
mode energy transfer rates (5^^(k'|p|q)) that represents the energy transfer rates from mode p of field X to mode k 
of field Y , with mode q acting as a mediator [irilll|. Note that k' -I- p -I- q = 0. This formalism is used for computing 
the shell-to-shell en ergy transfer because the earlier formalism suffers from ambiguity arising due to the third leg of 
the interaction (see The mode-to-mode energy transfer rates in MHD turbulence are given by 

^""(k'lplq) = -5([k'.u(q)][u(k')-u(p)]), 

5^^(k'|p|q) = -3([k'.u(q)][b(k')-b(p)]), 

^"''(k'lplq) = 3([k'.b(q)][u(k')-b(p)]), 

^''"(k'lplq) = 3([k'.b(q)][b(k')-u(p)]), 

where the above four formulas denote the energy transfers from u(p) to u(k), from b(p) to b(k), from b(p) to u(k), 
and from u(p) to b(k) respectively. For the derivation the reader is referred to the original papers (ifl [Tl]|. 

Using the above formulas we compute the shell-to-shell energy transfer in MHD turbulence by summing up the 
energy transfer among the Fourier modes. The energy transfer rates from m-th shell of field X {u or b) to n-th shell 
of field Y {u or b) is 

k'eri pGm 

The p-sum is over m-th shell, and the k'-sum is over n-th shell Since S'^"^(k'|p|q) is a fiuctuating quantity, we 

perform ensemble average to compute the average shell-to-shell energy transfer. We compute the ensemble average of 
S using a standard field-theoretic technique, a technique similar to EDQNM calculation fS, ES, fl^l . The calculation 
is quite standard, and it can found in McComb [l^l or Verma [Tll |. 

The function (S'^"^(k'|p|q)) depends on kinetic energy {E^{k)), magnetic energy {E''{k)), cross helicity {Hc{k)), 
magnetic helicity {HM{k)), kinetic helicity {Hxik)), and mean magnetic field. The total cross helicity, magnetic 
helicity, and kinetic helicity are defined as u-b/2, a • b/2, and u-[x)/2 respectively, where a and w are vector potential 
and vorticity respectively. The spectrum of these quantities are defined appropriately (refer to Verma for details) . 
For simplification we compute S for zero cross helicity and zero magnetic field, and in the inertial range using 
Kolmogorov's energy spectrum. To study the effects of kinetic and magnetic helicities, we have split S into helical and 



nonhelical components. The simplified expression (given below) is a function of Alfven ratio (va = (k) / (k)) , 
normalized magnetic helicity {vm = kHuik) / E^{k)), and normalized kinetic helicity [tk = HK{k)/ {kE^{k))). Please 
note that we are working in three dimensions. 
After lengthy algebra we obtain 
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[Pno^heUcal {v,w)+ F^J^^^i {v , w)] . (1) 



where JC" is Kolmogorov's constant for MHD turbulence, and some of the F^^j^^^^^i and F^J^-^^i are 
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Herei^*and 77* are renormalized viscosity and resistivity parameters. In this paper, we consider (a) nonhelical MHD 
{rM = tk — Q and different r-^s), and (b) helical MHD (r^ — I^tk = 0.1, tm — —0.1). For nonhelical MHD with 
TA = 0.5J,2^the constants {K'',v\ri*) taken are (0.55,2.1,0.5), (0.75,1.0,0.69), (1.0,0.64,0.77) respectively (see 
Verma (I6L l27l | for the procedure to compute these constants). For helical MHD, our choice of vk = 0.1 (small 
positive) and tm = —0.1 (small negative) is one of the typical values taken in numerical simulations, or observed in 
astrophysical situations; for this case, the constants if" = 0.78, v* = 1.0, ry* = 0.69 have been taken from Verma [2^ . 

The wavenumbers shells are binned logarithmically with the n-th shell being (/cqs""^ , fcos")- Note that the parameter 
s is similar to the scale-disparity parameter of Zhou We nondimensionalize the equations using the transformation 

m 

, a a a 

k = —; p = —v; q^—w, (6) 

u u u 

where a = fegs"^^. The resulting equation is 

YX 1 .u, /-i+t 

dw {VW) [F^onheUcal {v,w)+ F^J^^^i (w, w)] , (7) 
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which is independent of a in the inertial range. The independence of a implies self-similarity in the inertial range. 
From Eq. Q we can draw the following inferences: 

1. The shell-to-shell energy transfer rate is a function of n — m, that is, ^nm = ^{n-i){m-i)- Hence, the turbulent 
energy transfer rates in the inertial range are all self-similar. Note that this property holds only in the inertial 
range. 

2. T^^/n = —T^^/n, or b-to-u energy transfer rates from shell m to shell n is equal and opposite to the u-to6 
energy transfer rates from shell n to shell m. Hence T^^^/U can be obtained from T^j^/H by inversion at the 
origin. 

3. The MHD energy fluxes are related to the shell-to-shell energy transfers by the relationship 

n— m+l 

4. Net energy gained by a w-shell from u-to-u transfer is zero because of self similarity. However, a M-shell can gain 
or lose a net energy due to imbalance between u-to-5 and 5-to-u energy transfers. By deflnition, we can show 
that net energy gained by an inertial u-shell is 



{rpuh rpbu \ _|_ rpuh i q\ 
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Similarly, net energy gained by a 6-shell from b-to-b transfer is zero. However, net energy gained by an inertial 
6-shell due to u-to-b and b-to-u transfers is 
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Now we compute the integrals of Eq. Q; we denote the nonhelical and helical parts by (T^^f) ^ , , and 



nonhelical 

i'^nm) helical respectively. Their properties are described below. 



1. Nonhelical shell-to- shell energy transfer 

We compute the nonhelical shell-to-shell energy transfer using F^J^heiicai- We have chosen s = 2^/^^. This study 
has been done for various values of Alfven ratios. Fig. ^ contains plots of (TXJf /H) ... , vs. n — m for four 

^ ' — ' ^ V nm I } nonhelical 

typical values of ta = 0.5, 1, 5, 100. The numbers on the plots represent energy transfer rates from shell m to shells 
m-|-l, m + 2, ...in the right, and to shells m — l,m — 2,...in the left. For va = 0.5, the maxima of b-to-u energy transfers 
i'^nm) nonhelical ("^nm) nonhelical occurs at m = 71, and its values are approximately ±0.1 respectively. The 

corresponding values for r^i = 5 are approximately =f0.053. By observing the plots we find the following interesting 
patterns: 

1. The u-to-u energy transfer rate from shell m to shell n (T^m) nonhelical is positive for n > m, and negative for 
n < m. Hence, a u-shell gains energy from smaller wavenumber w-shells, and loses energy to higher wavenumber 
w-shells, implying that the energy cascade is forward. Also, the absolute maximum occurs for n = m± 1, hence 
the energy transfer is local. For kinetic energy dominated regime, s — 2^1"^ yields T^^JIl w 35%, similar to 
Kraichnan's Test Mean Field model (TFM) predictions Q. 

2. The b-to-b energy transfer rate T^^/Ii is positive for n > m, and negative for n < m, and maximum for 
n — m±l. Hence magnetic to magnetic energy transfer is forward and local. This result is consistent with the 
forward magnetic-to-magnetic cascade (H^^ > O) [Tllli^. 

3. For TA > 1 (kinetic energy dominated), kinetic to magnetic energy transfer rate (Tnm) nonhelical ^® positive 
most of the shells. For n — m< —30 or so, the value is small and negative. These transfers have been illustrated 
in Fig. El a). Using Eq. © we find that each u-shell loses a net kinetic energy to 5-shells, hence the turbulence 
is not steady. This phenomena is seen for all > 1 , and it could be one of the processes responsible for dynamo 
action. For s — 2^/'*, (t!'") ... , /H « 1.4, and the (T!!") ... , /II is positive for n > m — 1 and negative 

7 V fi-a / nonnelicaL ' ' V "-'"^^ nonhencat ' ^ — ^ 

otherwise. 

4. For TA = 0.5 (magnetically dominated), magnetic to kinetic energy transfer rate (Tnm) nonhelical positive 
for most of the shells (see Fig. QJ. For n — m < —30 or so, the value is small and negative. In addition, using 
Eq. © 

we find that each 6-shell loses a net magnetic energy to u-shells, hence the turbulence cannot be steady. 
This phenomena is seen for all va < 1- For s — 2^^^, {Tnm) nonhelical ^® positive for n > to — 1 and negative 
otherwise. 

5. The observations of (3) and (4) indicate that kinetic to magnetic or the reverse energy transfer rate almost 
vanishes near = 1. We believe that the evolution of MHD turbulence toward ta ~ ^ in both decaying and 
steady-state is due to the above reasons. For ^ 1, MHD turbulence is not steady. This result is similar to 
Pouquet et al.'s prediction of equipartition of kinetic and magnetic energy using EDQNM calculation An 
analogous result was discovered by Stribling and Matthaeus (2^ in the context of the Absolute Equilibrium 
Ensemble (AEE) theory. 

The steady-state value of in numerical simulations (e.g. DebHquy et al. Dar et al. Alexakis et al. 
0, Mininni et al. Q, and Haugen et al. [30]) and solar wind (e.g., Matthaeus and Goldstein ^1]) is around 
0.5-0.6. The difference is probably because the realistic fiows have more interactions than those discussed above, 
e.g., nonlocal coupling with forcing wavenumbers, coherent dissipative structures j^^l etc. 

6. When rA is not close to 1 (ta < 0.5 or rA > 5), u-to-b shell-to-shell transfer involves many neighbouring shells 
(see F ig, ni l- This observation implies that u—b energy transfer is somewhat nonlocal as predicted by Pouquet 
et al. [l^ . 

7. We compute energy fluxes using T,^, and flnd them to be the same as that computed by Verma ^HH^I- Hence 
both the results, flux and shell-to-shell energy transfer rates, are consistent with each other. 



Debliquy et al. 0| performed a decaying MHD turbulence simulations and computed the shell-to-shell energy transfer 
rates. The kinetic and magnetic heHcity was approximately zero. DebHquy et al.'s results show that the shell-to-shell 
energy transfers are forward and local. They also find that in the magnetically dominated MHD, the energy transfer 
from the same shell is from magnetic field to velocity field. Since the numerical simulations start with = 1, the 
energy transfer rates for kinetic-energy dominated regime is not known numerically. Our theoretical results are in 
general agreement with the numerical results of Debliquy et al. As an example, the numerical values of shell-to-shell 
energy transfer rates shown in Fig. 9 of DebHquy et al. @ {ta ~ 0.4) is similar to our theoretical results shown in 
Fig. nj A major difference between theoretical and numerical values are for T^''/Il where theoretical value for n = m 
is larger than its numerical counterpart. 

DebHquy et al. [6] found that T^" changes sign near ta ~ 0.5. In our theoretical calculation, the change of sign 
takes place around ~ 1. These findings are very encouraging, and they yield explanation why the asymptotic state 
of Alfven ratio is close to 0.5. The difference between theory and numerical simulations may be due to the neglect 
of large-scale non-Kolmogorov-like behaviour in theory, or due to the fact that Debliquy et al.'s results are based on 
decaying simulations, while the theoretical results assume steady-state. These issues need to be addressed. 

Alexakis et al. and Mininni et al. j3| computed shell-to-shell energy transfers in forced MHD. In the inertial 
range, the energy transfers are essentially local. However, the forcing velocity-shell (at large length-scale) provides 
energy to the large-scale magnetic field. The forcing velocity-shells also provide energy to the inertial range shells 
by nonlocal channel. Similar picture was observed by Dar et al. in forced 2D MHD turbulence. Our theoretical 
results are consistent with the Alexakis et al.'s results in the inertial range. Unfortunately, the present theoretical 
calculation cannot predict the coupling with the forcing shell. 

After the above discussion on nonheHcal MHD, we move to helical MHD. 

S. Helical Contributions 

Now we present computation of (Tnm) helical shell-to-sheH energy transfer rates for helical MHD [Hm 0, Hk 
0) fiSj. To simplify the equation, we consider only nonAlfvenic fiuctuations (cc = 0). We have chosen ta = l,rK = 
0.1, — —0.1. These values are one of the typical parameter values chosen in numerical simulations. We take 
s — 2^^^ to get a increased value for {^nm) helical ^'^^ above choice of parameters, Kolmogorov's constant 
— 0.78 f23|. In Fig. [2lwe have plotted {^nm) helical vs ti — m. Our results on helical shell-to-shell transfers are 
given below: 

1. Comparison of Fig. [21 with Fig. Q] (r^ — 1) shows that heHcal energy transfers are order-of-magnitude lower 
than the nonheHcal ones for the parameters chosen here (r^ = l,rK = 0.1, tm = —0.1). For maximal helicity, 
the helical and nonheHcal values become comparable. 

2. All the heHcal contributions are negative for n > m, and positive for n < m. Hence, helical transfers are from 
larger wavenumbers to smaller wavenumbers. This is consistent with the inverse cascade of energy due to helical 
contributions, as discussed by Pouquet et al. 

3. We find that the heHcal shell-to-shell energy transfer rate {Tnm) helical ^'^'^ i'^nm) helical significantly positive 
for —50 < n — m < 0. This signals a nonlocal b-to-b and u-to-b inverse energy transfers. Hence, helicity induces 
nonlocal e nerg y transfer between b-to-b and u-to-b wavenumber shells. This is in agreement with Pouquet et 
al.'s result [iJl that "residual heHcity" induces growth of large-scale magnetic field by nonlocal interactions. 

The theoretical findings Hsted above are consistent with Pouquet 's results based on EDQNM approximation and 
fiux calculations of Brandenburg et al. Alexakis et al. and Mininni et al. have computed shell-to-shell 

energy transfer in helical MHD turbulence; the helicity does change the energy transfer rates, however, in the absence 
of numerical value of normalized kinetic and magnetic helicity (r^c, r^/), we are not able to compare our results with 
their numerical values. 

In the next section we use our theoretical results and Debliquy et al.'s Q numerical results to argue why the 
asymptotic state of MHD fiows have « 0.4 — 0.6. 

III. CONNECTION WITH MHD ASYMPTOTIC STATE (ta ~ 0.4 - 0.6) 

The solar wind observations and numerical simulations show that the asymptotic state in the MHD fiows have 
TA « 0.4 — 0.6. We find in our theoretical analysis that for ta > 1, there is a preferential transfer of kinetic energy to 
magnetic energy; in fact, a u-shell loses a net amount of kinetic energy to 6-shell. For ta < 1, the pattern of u-to-b 



energy transfer is reversed, and there is a a net transfer of energy from magnetic to kinetic. This preferential energy 
transfer is minimum for ta ~ 1- Our theoretical calculation however is based on an assumption of Kolmogorov's 
spectrum for the energy, which is not valid for the smaller wavenumber modes. Here we use Debliquy et al.'s 
decaying simulations results for obtaining further insights into the large-scale energy transfers. 

Debliquy et al.'s f3| showed that the energy flux from small- wavenumber w— sphere to small wavenumber 6-sphere 
is positive for ta > 0.63, and becomes negative for lower r^. The global (inclusive of all shells) u-to-b energy transfer 
changes sign at r = 0.4. Hence, the smaller wavenumbers shells also play an important role in energy transfers. 
As a result, the asymptotic state of MHD turbulence has Alfven ration close to 0.5. Thus our detailed analysis of 
shell-to-shell energy transfer, and Debliquy et al.'s global and flux analysis is able to explain qualitatively why MHD 
turbulence evolves to asymptotic state with Alfven ratio ta ~ 0.4 — 0.6. Our theoretical predictions are consistent 
with the numerical results of Dar et al. whose asymptotic Alfven ratio ta is approximately 0.5. Alexakis et al. 
and Mininni et al. j3| also flnd their asymptotic Alfven ratio ta to be less than 1. 

IV. CONCLUSIONS 

In this paper we have computed the shell-to-shell energy transfers in MHD turbulence analytically. Our results 
provide theoretical explanation for the recently computed shell-to-shell energy transfers using direct numerical simu- 
lations. The contributions of nonhelical and helical terms have been calculated separately. We find that the nonhelical 
u-to-u and b-to-b shell-to-shell energy transfers are local as in fluid turbulence, i. e., most of energy from a wavenum- 
ber shell is transferred to the neighbouring shells. Comparatively, helical u-to-b and b-to-u energy transfers involves 
distant shells (nonlocal). 

We find that the helical shell-to-shell energy transfer is backward, that is from larger wavenumbers to smaller 
wavenumbers. For tk = 0.1, rjif = —0.1, one of the typical values observed in numerical simulations, the helical shell- 
to-shell energy transfer is order-of-magnitude smaller than the nonhelical ones. However, for maximal helicity, the 
helical shell-to-shell energy transfer is comparable to the nonhelical ones. In the present calculation, the parameters 
rx and vm have been chosen to be constants. This is a gross assumption considering that magnetic helicity and 
kinetic helicity have different signs at different scales. Even then we obtain results which are consistent with recent 
numerical results and earlier theories (Frisch et al. [131) • Hence, the present calculation appears to capture some of 
the essential features of helical MHD turbulence. 

Our results show that the inertial-range shell-to-shell for > 1, there is a preferential transfer of kinetic energy 
to magnetic energy; the direction of energy transfer switches for < 1. Debliquy et al.'s numerical simulations 
provide us important clues for the energy transfers at smaller wavenumbers. Using these results we can argue why 
the asymptotic state of MHD turbulence evolves to Alfven ratio of 0.5. 

Our theoretical results are in general agreement with the numerical results of Debliquy et al. Dar et al. fl(\ . 
Alexakis et al. Q and Mininni et al. 0, who observe local energy transfer in the inertial range. In forced MHD 
turbulence with large-scale velocity forcing, a significant energy transfer to large-scale magnetic field and nonlocal 
energy transfer to inertial-range magnetic field were found. Since the nonlinear energy transfers involves only the 
velocity and magnetic field variables, we expect that the features of energy exchange in the inertial range should 
remain approximately the same in decaying and forced MHD. However, the forcing at large-scale velocity shell would 
affect the large-scale magnetic field, and could also induce nonlocal interations. A theoretical model with forcing will 
be useful to understand forced MHD turbulence. 

Detailed pictures of energy transfer studied here is useful for understanding various physical process, for example, 
dynamo mechanism to generate magnetic field in astrophysical objects. These results are also useful in modelling 
MHD fiows and in simulations. For example, we need to model backscatter and forward energy transfer for large-eddy 
simulations. Refer to Debliquy et al. for discussion on backscatter and forward energy transfer in Fourier as well 
as real space. 

In EDQNM analysis, the wavenumber shells are logarithmically binned. In the present analysis of MHD turbulence, 
and many papers on fiuid turbulence show local energy transfers among wavenumber shells for nonhelical MHD. Hence, 
the local energy transfer assumptions made in EDQNM analysis is valid at least for nonlocal MHD j^]. The energy 
transfers are somewhat more complex for helical MHD. 

Pouquet et al. |Ti| performed extensive EDQNM analysis of MHD turbulence and showed that both local and 
nonlocal interactions exist in MHD turbulence. The local interactions cause the energy cascade, while the nonlocal 
ones cause equipartition of kinetic and magnetic energies. Pouquet et al. argued for inverse cascade of magnetic energy 
from the competition between helicity and Alfven effect. Our analytic calculation also predicts inverse magnetic-energy 
cascade due to helicity. Our calculation shows that near equipartition of magnetic and kinetic energy is due to complex 
process involving inertial-range shell-to-shell interactions and small wavenumber shells. This picture is some what 
different than that of Pouquet et al. 



To conclude, the shell-to-shell energy transfer rates provide important insights into inertial-range energy exchange 
processes in MHD turbulence. 
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Figure Captions 

Figure 1: The plots of shell-to-shell energy transfers {Tnm) nonhelical ^ ~ ^ for zero helicities (ctc = 
rR = rM = 0) and Alfven ratios = 0.5,1,4,100. Here s — 2^/'^'^. The u-to-u, b-to-u, u-to-b, and b-to-b are 
represented by dotted, dashed, chained, and solid lines respectively. For ta — 0.5, the maxima 

of {Tn^n)nonheUcal 

and (Tnm) nonhelical ±0.1 respectively, out of scale of the plot. The corresponding values for ta = 5 are =f0-053. 



Figure 2:Helical contributions to shell-to-shell energy transfers(r,^)^^j.^^j /II vs. n — m in helical MHD with 
TA = l,rK = 0.1, TM = -0.1 and = 0. Here s = 2^^'^. 



Figure 3: Schematic illustration of nonhelical shell-to-shell energy transfersT,|^y^/n in the inertial range for (a) 
kinetic-energy dominated regime, and (b) magnetic-energy dominated regime. In (a) the u-to-6 energy transfer rate 
T^/Ii is positive for n > m— 1, and negative otherwise, while in (b) the h-io-u energy transfer rate T^^/U is positive 
for n > m — 1, and negative otherwise. The u-to-u energy transfer rate T^^, and the 6-to-6 energy transfer rate 
are forward and local. 
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Figure 1: The plots of shell-to-shell energy transfers {^nm ) ^^^f^^^^^j /n vs. n — m for zero helicities (ctc = tk = tm = 0) 
and Alfven ratios ta = 0.5,1,4,100. Here s = 2^/=*^ The «-to-ti, &-to-M, it-to-b, and 6-to-ti are represented by dotted, 
dashed, chained, and solid lines respectively. For = 0.5, the maxima of [T-nm) ^^^f^^^^^i /H and (^^nm) /n are ±0.1 

respectively, out of scale of the plot. The corresponding values for rA~5 are =f0.053. 
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Figure 2: Helical contributions to shell-to-shell energy transfers(r,J^) ^^^^^^^ /H vs. n — m in helical MHD with = l,rK = 
0.1, rM = -0.1 and fXc = 0. Here s = 2^^*. 




Figure 3: Schematic illustration of nonhelical shcU-to-shell energy transfersT,^/n in the inertial range for (a) kinetic-energy 
dominated regime, and (b) magnetic-energy dominated regime. In (a) the M-to-6 energy transfer rate T^^/II is positive for 
most n, while in (b) the 6-to-w energy transfer rate T^^/II is positive for most n. The w-to-w energy transfer rate T^m, s-nd 
the 6-to-6 energy transfer rate T^m are forward and local. 



